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ABSTRACT 
We consider the multiplicity of stellar systems with (combined) magnitude brighter than 6.00 
in Hipparcos magnitudes. We identify 4559 such bright systems (including the Sun), and the 
frequencies of multiplicities 1,2,..., 7 are found to be 2718, 1437, 285, 86, 20, 11 and 2. We 
discuss the uncertainties, which are substantial. We also consider the distributions of periods 
of orbits and suborbits. We note that for even more restricted set of 478 systems with Ун « 
4.00, the proportions of higher multiples up to sextuple are progressively larger (213, 179, 
54, 19, 8, 5), suggesting substantial incompleteness in even the reasonably well studied larger 
sample. 

This sample can be seen as relatively thoroughly studied for multiplicity, and reasonably 
representative of stars more massive than the Sun. But the restriction to Ун < 6 means that 
our sample contains hardly any systems where all components are low-mass main-sequence 


stars (K or M). 


Data on multiplicity are important as a constraint on (1) the star formation problem, (ii) the 
problem of the evolution of the Galactic stellar population and (iii) the interaction of dynamics 
and evolution through the effect of Kozai cycles. We discuss these topics briefly. 


Key words: binaries: close — stars: statistics. 


1 INTRODUCTION 


The statistics of stellar multiplicity, i.e. the number of components, 
the distribution of periods, mass ratios etc., are poorly known, es- 
pecially for multiplicities higher than two. Yet it is important in 
many respects, for instance as a characteristic of the star formation 
process, and as an initial condition for stellar evolution. Our goal 
here is to determine the observed multiplicity of a reasonably well 
defined, well observed and moderately large bright star sample. This 
is probably the best-studied sample of its size, of stars sufficiently 
massive to be important for Galactic evolution. The disadvantages, 
as well as advantages, of using a magnitude-limited sample rather 
than distance-limited sample are discussed below. In a forthcoming 
paper, the observed statistics of multiple stars will be compared 
with a simulated sample, including selection effects. 

The bright star catalogue (BSC: Hoffleit & Jaschek 1983) is 
a fundamental resource when considering the stellar population 
of the Galaxy, or at any rate the nearer parts of the Galaxy. It 
lists multiplicities, but these are often visual multiplicities that 
may include line of sight, or ‘optical’, coincidences. Although 
roughly limited to magnitude 6.5, it is not entirely complete to 
this magnitude, and also includes several fainter stars. The multi- 
ple star catalogue (MSC: Tokovinin 1997) carefully identifies many 
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multiple systems, but restricts itself to multiplicity >3. The Hippar- 
cos catalogue (HIP: Perryman et al. 1997) has useful data such as 
parallaxes and proper motions that can help to distinguish optical 
from physical systems. The MSC provides Johnson V magnitudes 
of the brightest ‘resolved’ companion, rather than combined mag- 
nitudes for the whole systems. The MSC is constantly updated, and 
at the time of writing it contains 348 systems with multiplicity >3 
and У < 6. 

Before counting bright multiple systems, it is necessary to define 
both ‘system’ and ‘bright’. At a first glance, one would like to define 
a system as a collection of stars that are gravitationally bound to 
each other, but not to neighbouring systems. Unfortunately, because 
gravity is a long-range force, it is difficult if not impossible to draw a 
clear boundary. The entire Galaxy can be viewed as a single system. 
Intermediate between the Galactic scale and the scale of individual 
stars are Galactic clusters, globular clusters (at least two of which 
are bright enough in total to qualify) and various collections of stars, 
such as groups and associations, which might qualify as 'systems'. 
We discuss this issue in Section 2. 

For ‘bright’, we choose the HIP magnitude scale, as being rea- 
sonably homogeneous. But the issue of magnitude is somewhat 
complicated by the fact that we wish to use the combined magni- 
tude of the system. The 348 systems of the MSC referred to above 
involve Johnson rather than Hipparcos magnitudes, and when ad- 
justed for Hipparcos magnitudes the number is 330. 

Three comparable seventh magnitude components can make a 
sixth magnitude system. The obvious alternative, that we might use 
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the magnitude of the brightest component, seems to us to introduce 
an unnecessary extra uncertainty, since for any observing technique 
there will be systems that are marginally resolved, and the individual 
magnitudes will be less certain than the combined one. However, 
logic then dictates that we combine magnitudes even for systems 
that are well resolved. Some bright systems have components sev- 
eral hundred seconds apart. A part of our purpose is to compare 
the observed distribution of multiples with a theoretical model, and 
in the latter it is obviously most rational simply to combine the 
magnitudes of all the components. 

Our main reason for preferring the Hipparcos magnitude scale 
is that it averages the magnitude of variable stars in a logical and 
systematic way. The larger amplitude variables, such as ё Cep and 
Mira variables, are often listed in catalogues under a magnitude that 
is not a systematic average, and which may differ from the average 
by half a magnitude or more. 

Our main aim in this paper (Section 3) is to determine as well as 
possible the observed distribution of multiplicity, along with the dis- 
tribution of periods, for a reasonably large and yet reasonably well 
studied and complete sample. The fact that the sample is magnitude- 
limited rather than distance-limited makes it unrepresentative of the 
lower end of the mass spectrum, but is we believe compensated 
by the fact that there is in effect much greater signal-to-noise ratio 
(S/N). Somewhat by coincidence, this sample is in practice rather 
well representative of those stars massive enough to have significant 
evolution in the course of a Hubble time. A distance-limited sample 
would have to go out to ~250 pc to include any О stars; it would 
then include over 107 stars, the vast majority of which would not 
be at all well studied as to multiplicity. In addition, a considerable 
majority would be too low mass to have significant evolution within 
a Hubble time. 


2 MULTIPLICITY 


The great majority of multiple systems are ‘hierarchical’, with say 
a wide ‘binary’ containing two closer ‘binaries’. А few hierarchies 
can be as many as four deep, as in v Sco (HR 6026/7) which is 
one of the two septuple systems in the bright sample. Since there 
is usually a factor of ~107—10 difference in separation, from one 
level of hierarchy to the next, such four-deep systems are likely to 
be rare. This large factor accounts for the stability of hierarchical 
multiples, although in principle, and if the orbits and suborbits are 
roughly coplanar, a factor of only 3 or 4 allows long-term stability. 
However, non-hierarchical systems exist, in small numbers. They 
can be expected to be young and to break up in a few million years 
as a result of gravitational interaction, and presumably that is why 
they are rare. 

But there is a substantial grey area where non-hierarchical mul- 
tiples can overlap with clusters. Perhaps the entire Pleiades cluster 
should be seen as one non-hierarchical system. We prefer to see it 
as several independent hierarchical systems, but it is not clear that 
there is any sensible criterion that would compel us to do this. 

Note that we use the word ‘system’, and even the words ‘multiple 
system’, to include the possibility of multiplicity one, i.e. single 
stars. On the whole, we will avoid the word ‘star’ because this is 
often used ambiguously, to mean either a system if the individual 
components are not readily distinguishable, or a component of a 
system if they are. In this paper, we count systems, and within 
systems we count components, so far as we are able. We should add 
that we do not include planets as components, although we note 
systems that contain planets. This, of course, introduces another 
grey area, since it is not clear where stars stop and planets begin. 


The prototype of a non-hierarchical system is the Trapezium, 6! 
Ori (HR 1893-6), which consists (somewhat surprisingly) of five 
bright components arranged roughly on the boundary of an ellipse 
with axes 23 x 13 arcsec’: see the adaptive-optics mosaic of Simon, 
Close & Beck (1999, their fig. 1). The components are distributed 
non-uniformly around this ring, with angular separations (as seen 
from the centre of the ellipse) of about 90°, 30°, 30°, 120° and 90°, 
starting from the brightest component at the South and proceeding 
anticlockwise. But perhaps we should add to these five components 
the two main components of 0? Ori (HR 1897), about 50 arcsec apart 
and about 150 arcsec from the ellipse towards the Southeast. Further, 
these seven components have close subcomponents (12, according 
to Preibisch et al. 1999). One component, BM Ori (HR 1894), is 
in fact in a non-hierarchical quintuple subsystem: BM Ori itself, 
which is an eclipsing and spectroscopic binary, and two further 
components making a roughly equilateral triangle about 1 arcsec on 
the side, one of which is a very close pair separated by ~0.1 arcsec. 

However, we ought not to ignore several thousand other stars 
which are heavily concentrated towards the Trapezium, and evi- 
dently physically associated at some level. Several are within the 
elliptical curve (at least in projection), and not just outside it. We 
choose to consider the Orion Nebula Cluster (ONC) as just two 
independent bright systems, 0! Ori C (HR 1895) and 6? Ori A 
(HR 1897), which are the only two among the seven major com- 
ponents which qualify with Ун < 6. They are both hierarchical 
systems of multiplicity three, if we ignore their possible relation- 
ship to moderately close companions. 

We identify 17 non-hierarchical systems in our sample. Several 
of these may be simply the three or more brightest stars in a cluster; 
and there are a further few possible binaries that may actually be just 
the two brightest stars in a rather distant cluster (~ І kpc). Some hi- 
erarchical systems appear as non-hierarchical configurations simply 
because of projection. 

Among the bright star sample are several looser groupings, such 
as ‘OB associations’ and ‘moving groups’. These appear to be 
the last stages in the evaporation of stars from dense star-forming 
regions (SFRs) into the general Galactic field. They are probably 
not bound, and so we treat the members of these groupings as 
individual systems, but possibly multiple on a small scale, and if so 
then usually hierarchical. 


3 PROCEDURE 


We searched a number of catalogues, of which we have already 
mentioned three (the BSC, MSC and HIP). The BSC is partic- 
ularly helpful in view of its copious notes. We also used the 
eighth spectroscopic-binary catalogue (BFM; Batten, Fletcher & 
McCarthy 1989) and the ninth (SB9; Pourbaix et al. 2004). These 
two catalogues are also very helpful in view of their copious notes. 
We added the Catalogue of the Components of Double and Multiple 
Stars (CCDM; Dommanget & Nys 2002) of close multiples (many 
of which, however, are optical rather than physical multiples) and 
the sixth catalogue (VB6) of visual-binary orbits (Hartkopf et al. 
2000). We also incorporated information from the GCVS (Samus 
et al. 2004) on eclipsers and ellipsoidal variables, and from the 
CHARA survey of speckle observations (Hartkopf, McAlister & 
Franz 1989; Hartkopf et al. 2000). 

The Hipparcos experiment also generated a catalogue of dou- 
ble and multiple star measurements, HDMC. It contains mostly the 
known visual binaries, but also new systems discovered by Hippar- 
cos. HDMC frequently obtained solutions for pairs that gave the 
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same parallax and proper motions for close pairs, thus normally 
confirming their physical association. 

We included two tables of data on astrometric binaries by 
Makarov & Kaplan (2005, hereafter MK), where they considered 
systems whose Tycho and Hipparcos parallaxes and proper motions 
showed non-linear behaviour with time, i.e. acceleration, suggesting 
astrometric binaries. There are 348 systems containing such 'astro- 
metric accelerators’ in our catalogue. We also included tables from 
the Galactic Kinematics catalogue of Famaey et al. (2005, GKC ), 
who considered the space motions of a large number of HIP targets 
and noted those for which there was evidence of radial-velocity 
variations of an orbital character. Several have published orbits, but 
many more are from the private collection of R. F. Griffin. For sev- 
eral systems with composite spectra, consisting typically of a red 
giant and an A or late B star, we give the spectral types obtained 
by R. E. M. Griffin (soon to be published) as a result of careful dis- 
entangling of the two spectra; these types are often rather different 
from those normally quoted. 

In addition to GKC above, we used several catalogues of radial 
velocities (Andersen & Nordström 1983; Nordström & Andersen 
1983, 1985; Grenier et al. 1999a,b) for B, A/F and G/K/M stars in 
the southern sky. They used the results of several radial-velocity 
measurements per star to determine significant variation between 
observations. Although there were too few observations to establish 
orbits, they found several significant variables, identified as ‘VAR’. 
They also found several marginal cases, indicated as ‘VAR?’. We 
included the former as binaries (or subbinaries), but left out the 
latter. We included the catalogues of Dugennoy & Mayor (1991) 
who considered many multiple systems containing bright solar-type 
stars; of Harmanec (2001) for binaries including Be stars; of Aerts & 
Harmanec (2004) for binaries including pulsating stars; of Parsons 
(2004) for triple systems with cool giants and hot dwarfs; and of 
Lindroos (1985), who attempted to distinguish between physical 
and optical multiples among many bright and wide systems. 

Another catalogue was a private one maintained by one of us 
(PPE), with about 3000 entries taken from the literature in the 
interval 1975—2005; this catalogue concentrated on systems which 
contained stars evolved beyond (or not yet up to) the main sequence, 
but also contained many bright systems of main-sequence stars. We 
do not reference this catalogue directly, but instead refer directly 
to those papers in it, which supplied data that were different from 
(and, as we judged, better than) data from the principal catalogues 
mentioned above. 

Following McClure (1983) and Boffin & Jorissen (1988), we 
assume that all Ba stars are binaries, with a white dwarf (WD) 
component. Many Ba stars have indeed been determined to have 
spectroscopic orbits, and in a small number a WD has actually been 
detected in the ultraviolet. It is surprising to us, however, that among 
the 52 Ba stars in our bright sample only 10 have known orbits. 
Three of these have known WD components: £! Cet (HR 649), 
¢ Cyg (HR 8115) and ¢ Cap (HR 8204). However, the case for 
binarity is not just that some are confirmed binaries, but much more 
strongly that a physical mechanism exists to explain the Ba anomaly 
in terms of binarity, specifically with a WD companion, and that the 
anomaly is very hard to explain otherwise. 

The situation is rather different for S stars. Some of these may 
be evolved Ba stars, and so also with WD companions, whether 
seen or not. But others could be ‘intrinsic’ S stars, having produced 
their own abundance anomalies internally. We find only four S stars 
in our bright sample (including one with spectrum MAIIIS); two 
of these have known WD companions. The other two have visual 
companions, which are however too far away to be likely WD 
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remnants of S process donors. Perhaps, the systems are triple, but 
since they might instead contain intrinsic S stars we consider them 
for the present to be merely binary. 

For pairs of stars that are fairly close but might be optical rather 
than physical, we define a propinquity parameter 


X = log p + 0.3V — 3.95, (1) 


where p is the separation in arcsec and V is the magnitude of the 
fainter component. We expect X < 0 if the component is sufficiently 
bright and near that only ~ 1 such coincidence is likely in 5000 cases. 
The model is based on the approximation that the number of stars 
brighter than V is 


logn ~ 3.7 + 0.6(У — 6.0), Q) 


and that they are randomly distributed over the sky. Many pairs 
satisfy the propinquity test, and most of these are already fairly well 
established as either having measured visual orbits, or common 
proper motions. However, a handful is equivocal: they may satisfy 
the propinquity test by a considerable margin, and yet the proper 
motions from HIP or HDMC, and even the parallaxes, may be quite 
widely different. Often, though different, they are not significantly 
different because the measurement errors happen to be several hun- 
dred times larger than is normal; it is not clear why the errors are 
so much larger in a few cases. We tentatively identify a handful of 
systems where we suspect that the multiplicity is three rather than 
two, and that is responsible for (i) discrepant proper motions and 
(i1) unusually large errors. 

On the other hand, faint physical components with X > 1 are 
effectively lost in the stellar background, especially near the Galac- 
tic plane. Without further work on proper motions, e.g. Lépine & 
Bongiorno (2007), the propinquity test alone biases the observed 
statistics towards bright companions. 

We accept that some of the systems with negative X are optical, 
one (HR 2764) despite having X ^- —0.7. The parallaxes in this case 
are very different, but are roughly in agreement with the parallaxes 
expected from spectral types. The spectra are K3Ib-II and FOV, 
and yet the magnitudes differ by only 1.8. The parallaxes differ 
by a factor of 13, which is no doubt very uncertain since one is 
0.001 arcsec; but it seems more reasonable that the luminosities 
differ by a factor of 100-1000 than merely 5. A further optical pair 
(HR 6008/9) has X substantially less than zero; two exceptions in 
~5000 are acceptable statistically, just. 

Parsons (2004) has identified a substantial number of n > 3 
systems among stars whose spectral energy distribution (SED) has 
been measured over a wide wavelength range (0.13 — 0.9 u), putting 
together data from International Ultraviolet Explorer, Tycho and 
ground-based measures. The presence of two (at least) separate 
sources, one hot and one cool, is much more evident with such 
a range than in an SED that is limited to the classic UBV range 
(0.39 — 0.6 u). It is often possible to determine the two tempera- 
tures separately, by fitting the SED to judicious combinations of 
single-star SEDs, and this also determines the relative luminosities 
and radii. Since the parallaxes are known from Tycho/Hipparcos 
as well, the absolute luminosities and radii are known, and can be 
compared to theoretical isochrones. Parsons finds several systems 
where the hot star is too bright, compared with a theoretical model, 
by a factor of —2, and concludes that the hot component is a sub- 
binary of two roughly equal components. Parsons (2004) lists 19 
systems which feature in our catalogue. In four of them, the hot 
component is already known to be a subbinary of short period, ei- 
ther spectroscopically or photometrically; although in one of these 
four (В Cap, HR 7776), the subbinary is single lined and therefore 
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the unseen subcomponent cannot be contributing significantly to 
the SED. However, HR 7776 (Table 1) appears to be sextuple, and 
it is possible that the nearest of the three remaining components, at 
~0.8 arcsec, contributes something. 

In estimating the multiplicity for a particular system in the present 
compilation, we use a different, and weaker, criterion than the MSC 
for the certainty with which the multiplicity has been determined 
so far: in legal terms, our criterion is roughly ‘on the balance of 
probabilities’ rather than ‘beyond reasonable doubt’. Given our 
present knowledge, we ask what is the most probable multiplicity. 
Clearly in some marginal cases, the multiplicity with the highest 
probability may have a probability only slightly over 50 per cent. 
For inclusion in the MSC, the criterion was normally stronger; that 
radial-velocity orbits had actually been determined, for example, 
rather than suspected on the grounds of significant radial-velocity 
variability. 


4 RESULTS 


Table 1 illustrates our results. The main body of the table is avail- 
able electronically in full; and in the printed version here, Table 1 
includes only a few examples, with a range of multiplicities. Many 
systems consist of two or more HR entries; we list them under the 
largest relevant HR number (Column 1). Table 2 is a table of cross- 
referenced names (not shown here in full but available electronically 
as Supporting Information in the online version of this article); we 
give the HR, HD, HIP and other identifiers for all components that 
have HR numbers, so multiple systems which include several HR 
numbers can be easily found. For systems at the end, which have 
no HR number, we give a ‘pseudo-HR’ number greater than 9200 
and prefixed by P. These are identified with genuine names (HD, 
HIP, ...) in the same cross-reference table. 

Column 2 of Table 1 lists what we consider to be the most rea- 
sonable multiplicity. A multiplicity n that is 73, and not followed 
by a query or colon, is the same as in the MSC (316 cases); a colon 
indicates that it is in the MSC with a different n > 3 (13 cases) and a 
query indicates that it is not in the MSC (75 cases). Most of the last 
category are systems where we feel that the balance of probability 
lies with n 7 3 but that the information is not so compelling as to 
make the multiplicity say 95 per cent certain. 

Column 3 contains some reference letters, defined in more de- 
tail shortly. Column 4 contains the parallax from Hipparcos, except 
that a Hipparcos parallax «0.001, including values like —0.002, is 
replaced by 0.001 automatically. Column 5 contains our descrip- 
tion of the configuration of the system using nested parentheses in 
roughly the format suggested by Evans (1977). For each individual 
component, we give a magnitude and a spectral type, where we can 
find them, and for each pair of components we give either a period 
or a separation in arcsec, where we can find them. Our reason for 
preferring this notation here is that, from experience, it can sum- 
marize a system sufficiently clearly that one can readily see where 
each component is in the hierarchy, and yet it only takes one line 
per system. To convey substantially more information, both about 
the subcomponents and the suborbits, we would need at least one 
line for each subsystem, as in the MSC. 

Although the one-line summaries exemplified in Table 1 may not 
look machine readable, they are. A short code along with the data 
allowed the following Tables 4—6 to be generated automatically. 
Various subsets, e.g. systems with WDs, or with semidetached sub- 
binaries, can be quickly identified. 

The magnitudes come from HIP if available, and the spectral 
types generally from the BSC or the MSC if available. HIP magni- 


tudes are given to two decimal places (dp) and other magnitudes, 
mostly Johnson V, to only one. Sometimes combined magnitudes of 
subsystems are given as well, and always the overall combined mag- 
nitude. Sometimes we have to combine Johnson magnitudes with 
HIP magnitudes to reach a combined magnitude; but the Johnson 
magnitudes are usually quite faint and so we treat the combination 
as a HIP magnitude. 

The separation is given to 3 dp if it comes from Hipparcos, and 
only 2, 1 or 0 if from another source. The period is in days if it 
comes from a spectroscopic or eclipsing binary; the reference is to 
SB9 or GCVS by default, but the letter ‘s’ implies an individual 
reference listed in the electronic version. The period is in years if it 
comes from a visual binary; the reference is to VB6 by default. Most 
systems with multiplicity n > 3 come from the MSC, by default, 
but the letter ‘s’ implies an individual reference, and in some cases 
letters such as AE (see below) explain why we think it is triple, 
though not with as high probability as would warrant inclusion in 
the MSC. 

Some published visual orbits have periods in excess of 300 yr. 
We characterize all of these by a separation rather than a period, 
feeling that an orbit should be seen revolving at least once (and 
preferably twice) before being considered reliably determined. In 
fact, several visual orbits with periods substantially less that 300 yr 
are quite tentative. 

The reference letters in column 3 have the following significance. 
By default, i.e. in the absence of a reference letter, a system with 
n> 3comes from the MSC, and with п < 2 from SB9, USNO, BSC, 
GCVS and/or CCDM, supplemented by HIP or HDMC. Systems 
with orbits from SB9/GCVS or VB6 are distinguished by having 
the period in days or years, respectively. When a reference letter 15 
given it has the following meaning. 


(i) A is for an entry that has been identified as a probable astro- 
metric binary by MK; this presumably has a long period, and so 
if the system is already known to have a short period (from 589, 
GCVS, or another source), we assume the system is triple (at least). 
If this makes it a triple that is not in the MSC, we write n — 3? 
rather than n — 3. Quite often it is in the MSC, however, because 
there is convincing additional evidence. 

(ii) b is for an entry of two or more stars that are (usually) intrin- 
sically bright, distant and not very close together. They are close 
enough together that their juxtaposition is unlikely to be chance, yet 
far enough apart that one can question whether they are in a long- 
lived orbit. It ends up as largely a subjective matter whether one 
thinks these should be listed as one system, or two or more systems. 
In some cases, they may be members of a cluster, and (arguably) 
not much closer together than cluster stars usually are. We list eight 
systems that we qualify with b. 

(iii) C is for an entry that comes from CHARA. 

(iv) E is for an entry identified in the GCVS as having an eclips- 
ing, or eclipse related, light curve. If a system merits both A and E 
we regard it as a probable triple, since an orbit that gives eclipses 
will usually be too small to also give an astrometric indication. 
Most but not all systems flagged E from the GCVS have known 
periods, coinciding with 589 periods, but those without periods are 
presumably rather uncertain. 

(v) G is for a system whose parameters are taken from the cata- 
logue of R. E. M. Griffin (to be published). These are mostly G/K 
giants plus B/A dwarfs; but in a handful of cases, she concludes that 
the dwarf is itself in a subbinary. 

(vi) H is fora system that we identify, through slightly discrepant 
proper motions, as having an extra, unseen, component making 


No claim to original US government works. Journal compilation © 2008 RAS, MNRAS 389, 869-879 


ӯс0с AeW | у uo зѕәпб Ла 001 26/698/2/682/91о0е/зезишу/шоодпо'опшерезе//:зауц Wo. ререојимод 


873 


Bright multiple systems 


Downloaded from https://academic.oup.com/mnras/article/389/2/869/974100 by guest on 11 May 2024 


"зопалојол оф Aq ројеотрш se ‘uses Ароолр зои Inq ‘PONIJU ST jeu) JIGIO Ue 10 jueuoduroo v Jo зопозола оф sojeorpur (UOTRINSYUOD) с uumgo;) ш yeu uonsenb y 
(91/11/91 МН '8'2) sexeqpered pasy Jup 
јецмәшоѕ әлец 124 pue urojsÁs owes оф JO Wed әд ој Jwadde sjueuoduroo злош 10 OM) JI гразелоле (q) pue (зопјел әлце8әи pue олог зшрпјош) sty) uey, 5591 JI 2oso 1000 01 dn рорипол (е) Inq 'зооледфн 
wo xej[e1red оф sr р иштоо хәз оф ur pequosep ose se *ѕәолпоѕ опдојејео летпоплеа sodu ose IINA зополојол V JO oouesqe ou] хәз оф ш pequosep se '5аолпов SNONVA 0} JOJoI SION] оф *c UUNTOD up 
попа ујиш epqeqoud јѕош рәјешцѕә оф st (и) с ишпјод 'epdures Ino лој Ajpenb 

jeu 5шој5л8 YH-Opnəsd полов ATUO Ајпизрт әм “G9TTE рие 96706 dIH S! ‘€O7T6d ‘OY UMOYS ојашехо DUO ‘с әв], oouo19ja1-ssoJ5 оф ur punoj eq UI иво ѕләдшпи qH Jo/pue ан 8шриойѕәлшоо eur а 4q 
рохцола 1026 зодшти ,YH-Opnəsd, в олт8 ом *urojsÁs оф ur зәдшпи YH UMUS ou 1 олод у (І ошод) (мни, pere? зодшти мн [eurrxeur оф әѕп ом *jueuoduioo YH euo Uey, злош Surumejuoo ѕшә]565 104 
“QPL 510) JO UOISIOA оштио oY} ur попецџојиј SuTIOddng eu ur рипој әд џез әде) оф JO UOTSIOA PN} AJL 


IeATIISIA6' S 100'0 1 :L0c6d 

[ооволе Q6t7 6C I ‘AG/SAL9'9 + (i 2 + ЛОС 9] 6L S 700'0 ан ¿E :£0c6d 

МАСТ 1€0'0 E él 1106 

{әәѕолә g'coc (295018 89'0 1УРГ6 +ШОУ91'9)60%9 + [295022 8'0 ‘E'8 + (c0 = ә “Рӯи (960 = 2 *р8938 55 + ASMTL + IISD)IICE) РГЕ 600'0 оа 9 9144, 
(89'0 = 2 'PIOCC А-МОЗИ8 + ПЕМ) LLS с00'0 svu T :9709 

€'O-eJIIIS TAT T'O 680'0 5 1 ‘Ores 

(ооволе 1'67 *IIIOMS IT + 960) Les 200'0 T iT ‘8067 

[2?eso: 000 “100 Ф 68459 + ооволе LOT STE O (o ‘i + ошовор' + | “i + euATCEDCS Cl TET 8000 qvv € 1797 
[ооволе 00210 :ЛІ-Ш69#С/ + (o ‘i + A88D0OC'9] 16'5 L00°0 suu "s '£96€ 

[06'0 = ә 6886 4 + (AS "POET T 'A-AISD + ATH] 61:5 €c0'0 sav € ‘8817 

[ооволе с69'1 :ACH68 t + (CO = 2 05:6 (prog 0 : + :9) + (10:0 = ә 'Р066' 1, “OTH + АТ80)85 €l 6с 700'0 ЯЯ $ és 188/1 
[Sesame OOS TI :( + лӮ9)659 + (o + + ATOHL9'9)] 87'S 120'0 Suu t 9981 

(2:115 + судм) FL v 900'0 5 T :OSST 

(PO'STI ТИХ + IIO3D 89°¢ +00°0 sou С :856 

(295218 969:9] ‘UATVS6'9 + л69/0'9) LOS 210'0 Ч iT :679 

(68'0 = 2 'PO LOT ::00 + ASV) OL'T с60'0 5 c 158$ 

(әәоле [0 :л69 + 90У) 78'S 200'0 о (d гбЕР 

(ооволе [90 p£ T :919980`/ + #10411'9) 56% 100°0 q ic ‘TSE 

[£€'0 = 2 'PI60Z ‘(P9161 : + A69) + *IIISD)] LES +00°0 D ¿E єс 

[оевоје [/8c CIAO ЕТ + (Е'0 = 2 08510: + MEM] ЕРЕ «00 vu € :691 

(O€ 0 = 2 'р09/5 ‘i + ШОУ) 925 с00'0 u T 1081 

[9170 = 2 40689 :л0006'9 + (PC80'C ‘i + ASHI9'S] TES 8T0'0 € у 

{oasore opc :(»osoxe 1'0 ‘AOV + AOV)60'S + [998918 090' Lc ‘(PL0 = ә одру SALW + ACV)SS Y + (ооволе pc СЕТ + AGMEEH189'E} cre 100 9 :9ЄТ 
(i fi + ATA) SL'S 200 V iT ‘071 

[070 = 2 ALTSI ‘8'9 + (TTO = ? 'PCE' SC ‘i + AISS)S6 S] $55 200'0 5 € 16 

MISOIL’S 600'0 1 Y 

uoneinsyuod хе[елед WAPA u "Hu 


'suiojs&s ojdnjnur зу8џа Jo suoneun3guoo ојашес ‘T әде, 


No claim to original US government works. Journal compilation © 2008 RAS, MNRAS 389, 869-879 


874 


Table 2. Sample from cross-reference table. 


P. P. Eggleton and A. A. Tokovinin 


IDS/MSC Cluster 


CCDM 
/WDS 


ADS 


BD/CoD 


Flam- 


mHR HR HD HiP GCVS Bayer 


Galactic latitude 


Galactic longitude 


/CpD 


Steed 


00270-6331 


CP-63d50 
CP-63d50 
BD+62d102 


betTucl 


HP2484 


HD2884 
HD2885 
HD2905 


HR126 


HR136 


—54.02 
—54.01 
0.14 
—42.36 


—55.61 
—53.97 


—39.26 
32:82 


306.78 


00270-6331 


00315-6257 


НР2487 Бе! Тис2 


HR127 


HR136 


306.78 


Сазов14 


15Са5 
52Psc 
51Psc 


kapCas kapCas 


HP2599 


HR130 


HR130 


120.84 


ADS452 
ADS449 


BD+19d79 
BD+06d64 
CP-63d52 
BD-02d5631 


HP2568 
CHP2548 


HD2910 
HD2913 
HD3003 


HR131 


HR131 


116.94 


00272--0624 
00270-6331 


00324--0657 


NSV15113 


HR132 


HR132 


114.55 


betTuc3 


HP2578 
HP107516 


HR136 
PHR9207 


HR136 
PHR9207 


306.54 
54.20 


EPAqr 
NSVA171 


HD207076 


NGC2632 


08399+1933 
08399+1933 
08399+1933 


BD+20d2150 ADS6915 


BD+20d2152 
BD+20d2153 


PHR9209 HD73598 HP42497 


PHR9209 


206.01 


NGC2632 


ADS6915 


32.34 PHR9209 PHR9209  HD73618 


32.34 


206.00 


NGC2632 


ADS6915 


HD73619 NSV4174 


PHR9209 


PHR9209 


206.02 


The full version of this table can be found in the Supporting Information in the online version of this article. 


The first two columns are galactic latitude and longitude. The next two columns, mHR and HR, give the system number and the component number, which may differ if a system contains more than one component 


with an HR number – for example, HR 126/127/136. In a small number of cases at the end, systems and components have по HR number, and so are given a ‘pseudo-HR’ number, of no significance except to 
cross-reference Tables 1-3. The remaining columns are standard identifiers, except that Greek letters in Bayer names are rendered by the first three letters (or if necessary only two) of the usual spelling in English. 


the system at least triple. There are six such systems; we assign their 
multiplicities as 3?, indicating that we think that the probability of 
the extra body is over 50 per cent, though nothing like certain 
enough for inclusion in the MSC. 

(vii) h is for a binary system where we suspect that large but 
uncertain discrepancies in the HIP values of parallax and/or proper 
motion may be hinting very weakly at a third body. We assign 
multiplicity 2?, and we identify 25 such systems. 

(viii) Lis for a system where we have relied on the discussion by 
Lindroos (1985) as to whether certain companions, usually rather 
distant, are optical or physical, and arguably pre-main sequence if 
physical. 

(ix) P is for a system noted by Parsons (2004) as containing a hot 
component that may be two roughly equal components. 

(x) R is for an entry identified as a probable or certain spectro- 
scopic binary by Famaey et al. (2005), Andersen & Nordstróm 
(1983), Nordstróm & Andersen (1983, 1985) or Grenier et al. 
(1999a,b). 

(xi) s is for a specific reference which is given in the electronic 
version of the table. 


Although nested parentheses are a good way of displaying hier- 
archical character, they do not work for non-hierarchical systems. 
For these we use the following notation. Suppose that X, Y and Z 
are a non-hierarchical triple. Then, we write it as 


X + Y G6, р+2@0', р 


with position angle 0 in degrees and separation p in arcsec. Both 
are measured from component X. An example in the printed Table 1 
is HR 4621. 

We attach considerable weight to Hipparcos observations, but 
we should note that Hipparcos claims that а Boo (HR 5340), at a 
distance of only 11 pc, is a binary with AVy = 3.33 and separation 
0.260 arcsec. Griffin (1998) has shown that this is very difficult to 
reconcile with the record of radial-velocity measurements, which 
are many and accurate. It is hardly conceivable that such a compan- 
ion would not have produced a recognizable radial-velocity curve 
over the last 50 yr. Nevertheless, the issue may still be open: see 
Sóderhjelm & Mignard (1998) and Verhoelst et al. (2005). We list 
it as single, believing that the companion is an artefact. 

We discuss the small subset of entries listed in Table 1, mainly to 
aid the reader in seeing why we put forward the multiplicities listed 
there. We refer to the electronic version of the table as the ‘EV’. 


HR 4 The absence of reference letters implies that the parallax 
and magnitude come from HIP, and the spectral type from BSC; 
and there is no clear, or even tentative, indication of any companion 
in any of the catalogues consulted, or indeed in any paper we have 
seen. 

HR 91 An absence of reference letters, coupled with the fact that 
the multiplicity is n > 3, would mean that the basic source of the 
data is the MSC; but the magnitudes come from HIP. However, the 
reference letter ‘s’ means that some data (the shorter orbit) comes 
from a paper referenced in the EV, and in fact the longer period, 
taken from VB6, is moderately different from the MSC. 


HR 120 This is an ‘astrometric accelerator’ from MK, which we 
therefore suppose to be a binary with an unseen companion and 
unknown (but presumably fairly long period) orbit. 

HR 136 This is actually three HR stars (126/127/136), but we la- 
bel it by the largest HR number. The system is a wide but apparently 
hierarchical triple of three close binaries. The three HIP parallaxes 
differ by a surprising 20 per cent, and the three proper motions are 
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Table 3. Sample from reference table. 


Bright multiple systems 875 


mHR Reference 

91 ApJ655,473 

553 Obs 108,228 

958 JAA11,491 

965 Obs109,180 http://www.eso.org/gen-fac/meetings/ms2005/griffin.pdf 
1556 ApJ327,214 


The full version of the table can be found in the Supporting Information in the online version of this article. 


mHR is the maximal HR number as in Table 1. 


Table 4. Multiplicity frequency in four samples. 


Sample Total n=1 2 3 4 5 6 7 Average 
Vy <6 4559 2718 1437 285 86 20 1 2 1.53 
Ун<4 478 213 179 54 19 8 5 0 1.84 
North 2141 1233 697 140 52 11 7 1 1.57 
South 2417 1484 740 145 34 9 4 1 1.49 


The last column is the average multiplicity or ‘companion star fraction 
(CSF)’, defined as У), nN4/ У), Nn. If Nn c a^", the average is a/(a — 
1). The first two samples include the Sun, the last two exclude it. 


somewhat different, though not alarmingly so. Probably, the inter- 
nal motions of the three binaries account for most or all of this. The 
system is accepted as real by the MSC. Our parallax is the mean of 
the three. 

HR 142 With no reference letter (unlike HR 91), the data are 
straight from the MSC, except for the magnitude from Hipparcos. 

HR 152 This is a radial-velocity variable flagged by GKC; in 
fact it has a known orbit, in SB9 (by default), which we presume is 
responsible for the variation. 

HR 165 This is both a radial-velocity variable from GKC, and 
an astrometric accelerator from MK. We assume that both of these 
are accounted for by the spectroscopic orbit which is (by default) 
in SB9. The separation of the third component, quoted to only 1 dp, 
comes from the BSC. However, if the separation was quoted to 3 dp, 
it would be a HIP measurement, where there is normally supporting 
evidence from common, or nearly common, proper motion. The 
unseen companion in the spectroscopic orbit might, in principle, be 
a WD, but the mass function is sufficiently small that an M dwarf 
cannot be ruled out. 

HR 233 This n > 3 system is not in the MSC, hence '37. R. E. 
M. Griffin finds an inner orbit of 1.916 d. ‘G’ isa reference to R. E. 
M. Griffin (to be published). 

HR 382 A pair of intrinsically bright, distant stars, part of the 
cluster NGC 457. They are the two brightest members by ~2 mag, 
and some way to one side of the cluster centre as defined by about 
half-a-dozen stars of magnitude 9-10. Their projected separation 
is ~0.67 pc, about 2.5 pc from the centre. They seem to be on the 
rather broad margin between indisputable binaries, and pairs or 
small groups of stars that clearly have a common origin but may 
no longer be bound. We note eight such systems, referenced with 
letter b. 

HR 439 The separation comes from CHARA. Since it is small, it 
may be well changing quite rapidly; we simply list the value quoted 
at the epoch of measurement. 

HR 553 Although the orbit is in 589, the spectral type of the 
companion is not; the ‘s’ means that a specific reference is given 
for this. 

HR 629 The 'h' implies that this is a system where the difference 
in HIP proper motions is rather larger than we would expect as a 


result of the well-known companion at 17 arcsec. However, one of 
the proper motions is unusually uncertain. We therefore entertain 
the possibility that a third body is present to disturb one of the two 
known bodies, but do not feel that the probability is high enough to 
raise the multiplicity from 2? to 3?. 

HR 958 This is a radial-velocity variable from the GKC, but an 
orbit is known which presumably accounts for this. It is not in the 
589 so the flag ‘s’ implies a specific reference to it (in the EV). 
The spectral types are referenced to 'G' (i.e. R. E. M. Griffin, to be 
published); they are more precise than the K1IIIep + A6V given in 
the BSC. 

HR 1556 The ‘s’ implies a specific reference (in the EV), to the 
fact that the companion is known to be a WD. 

HR 1564 This visual binary contains two radial-velocity vari- 
ables, as indicated by the MSC (as well as by the double entry, 
implied by RR, in the GKC). 

HR 1788 This complex system is based on a specific reference 
(De May et al. 1996), as implied by ‘s’. In addition to the eclipse 
period of the first subbinary listed, there is a second period (0.864 d) 
which comes from some other component. It is most probably ellip- 
soidal variation, and while it might come from the B2V component, 
the second-last component listed, De May et al. conclude that it is 
probably from the component we have indicated. 

HR 2788 This is a well-known Algol (R СМа), but it is also an 
astrometric accelerator (MK). The latter can hardly be due to the 
small orbit of the Algol, and so we could infer a distant third body. 
In fact such a body is known, and listed. There are other systems 
where both ‘A’ and 'Е' make it clear why we suppose them to be 
triple. There are four such AE triples in our catalogue. SD in the 
description stands for ‘semidetached’; we note 18 such systems in 
the entire sample. We flag contact systems with CT, and there are 
three such systems. 

HR 3963 Hipparcos gives very different parallaxes (0.007, 
0.001 arcsec) and fairly different proper motions; thus, the system 
appears to be optical. Yet a factor of 7 in distance seems entirely at 
odds with the fact that the spectral types and magnitudes are fairly 
consistent with their being at the same distance. Also they are very 
close for an optical pair; we would only expect one such pair in 
a sample 10 times larger. We suggest that the smaller parallax is 
wrong because of an unseen third body, and that this accounts for 
the difference in proper motion also. 


HR 4621 This is a non-hierarchical triple in the MSC, but two 
components are astrometric accelerators. Thus, we conclude that it 
may have n — 5; the triple is too wide to account for the acceleration 
without other bodies. 

HR 4908 This apparently improbable juxtaposition of an O su- 
pergiant with a faint K giant companion is suggested by Lindroos 
(1985) to be a physical system in which the secondary is still con- 
tracting to the ZAMS. Lindroos suggests several other pre-Main- 
Sequence companions to young massive stars. 
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876 P.P. Eggleton and A. А. Tokovinin 

HR 5340 We have already discussed the fact that this star 
(а Boo) is perceived as binary by HIP, but can hardly be binary 
in the face of much radial-velocity information to the contrary. The 
‘s’ flag indicates that we give specific references. 

HR 6046 This is both an astrometric and a radial-velocity vari- 
able, but both are presumably accounted for by the known orbit. 
SB9 gives an approximation to the orbit, but our flag ‘s’ points to a 
specific reference (in the EV) which we believe is better. 

HR 7776 This sextuple is rather complex, with conflicting in- 
terpretations. We accept that of the MSC. The flag ‘P’ indicates 
that Parsons (2004) considered that the system needs an extra hot 
component. 

HR 9072 The flag ‘s’ points to two specific references which 
flatly contradict each other. One says the system is triple, the other 
single (or at least not demonstrably multiple). It is probably a y 
Dor variable, with intrinsic pulsation capable of masking, or being 
mistaken for orbital motion. 

P9203 This is a pair of non-HR stars (HIP 32256 and 32269) 
which are rather far apart on the sky and yet which are both suf- 
ficiently bright that the propinquity factor (X = 0.2) is uncomfort- 
ably small for an optical pair. Their parallaxes are the same, but 
their proper motions differ modestly, to the extent that we invoke an 
unseen third body to account for the discrepancy. Their luminosi- 
ties are uncomfortably equal considering the difference in spectral 
types, but not unreasonably equal. If they are a real pair, their com- 
bined magnitude puts them in our catalogue, and since they would 
then be an п > 3 system we flag them with ‘n = 37, i.e. not in 


Table 5. Period distribution in systems and subsystems. 


the MSC. We assign а ‘pseudo-HR’ number (>9201), and give the 
actual ID in the electronic cross-reference table. 

P9207 This is a Mira variable whose average magnitude as de- 
termined by HIP puts it in our sample. It is not in the HR catalogue, 
but is HIP 107516 as given in the cross-reference table. 


Table 4 gives the numbers of systems observed with multiplicity 
from one to seven, for four samples. The first two are Vy < 6 and 4; 
the last two are the first sample divided into northern and southern 
hemispheres. Note that the Sun is included as a bright star in the 
first two samples, but is not included in either the north or the south 
sample. 

Table 5 gives the distribution of orbital periods, in bins of width 
1.0 in log P (yr). For wide systems with only a separation listed, 
the period is estimated from the angular separation, parallax and 
Kepler’s law, assuming a mass dependent on spectral type. The 
distribution is broken up into separate distributions according to the 
spectral type of the main component; type B is divided into early 
B (BO to B3.5) and late B. The distribution of period can be seen 
to be strongly bimodal at early types, becoming weakly bimodal or 
unimodal with a very broad maximum, at later types. The minimum 
where it appears is generally in the bins 0.1—10 yr. 

Table 6 gives the fractional multiplicity, broken down by the 
spectral type of the leading component. We see that O stars (which 
for our purposes include two Wolf-Rayet systems) appear to be 
substantially more multiple than later types. However, there seems 
to be little variation in fractional multiplicity among types B-G. 


log P (yr) 3.0 2.0 1.0 0.0 1.0 2.0 3.0 4.0 5.0 6.0 7.0 8.0 Total 
spectral 

о 0 5 11 5 0 4 3 6 12 1 0 0 0 47 
eB 0 25 41 21 14 21 26 29 33 21 3 0 1 235 
IB 0 25 53 24 20 43 62 63 48 16 8 0 0 362 
A 0 27 62 25 46 78 78 66 47 24 9 0 0 462 
F 0 14 33 24 39 46 61 44 26 14 3 1 0 305 
G 1 7 9 20 40 38 49 38 32 23 6 1 0 264 
K 0 4 4 12 56 35 35 42 35 26 4 0 0 253 
M 1 0 0 3 7 6 4 9 9 3 3 0 0 45 
Sum 2 107 213 134 222 271 318 297 242 128 36 2 1 1973 


The first column gives the spectral type of the dominant body in the system: eB means early B, i.e. ВО – B3.5, and IB means later B. Wolf-Rayet 


stars (2) are included under O; S and C stars are included under M. The first column of integers gives the number of systems and subsystems 
with log P (yr) < —3.0; the second for —3.0 « log P « —2.0, etc. The two shortest periods are a contact binary subsystem of the О dwarf 44 
Boo (HR 5618), and a cataclysmic binary subsystem of the M giant CQ Dra (HR 4765). Long periods are estimates from the angular separation, 


distance and Kepler's law. 


Table 6. Fractional multiplicity by spectral type. 


Spectral 1 2 3 4 5 6 7 Average totg tot4 
[9] 0.342 0.263 0.184 0.105 0.053 0.053 0.000 2.42 38 7 

eB 0.531 0.312 0.115 0.020 0.015 0.002 0.005 1.70 401 84 
IB 0.541 0.329 0.093 0.032 0.002 0.003 0.000 1.64 653 63 
A 0.540 0.362 0.068 0.018 0.006 0.005 0.000 1.61 928 90 
F 0.526 0.379 0.062 0.026 0.007 0.000 0.000 1.61 578 49 
G 0.550 0.370 0.056 0.020 0.002 0.002 0.000 1.56 588 60 
K 0.709 0.253 0.031 0.008 0.000 0.000 0.000 1.34 1043 99 
M 0.821 0.155 0.021 0.003 0.000 0.000 0.000 1.21 330 26 


Early B stars (B0—B3.5) are called ‘eB’; later B stars are called ‘IB’. Wolf-Rayet stars are included under О; S and C stars under M. The third 
last column is the average multiplicity; the last two columns are the total number, in the larger sample (Ун < 6) and the smaller sample (Ун < 


4), respectively. 
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It is not surprising that our G/K sample should be little different 
from our A/F sample, since most of our G/K sample are evolved 
giants and are essentially the same population as the A/F sample 
(mostly main sequence); but there is little variation even for eB and 
1B, where the relation is less close. 


5 COMPLETENESS, DETECTION EFFICIENCY 
AND SELECTION EFFECTS 


One might suppose that it would not be difficult to identify the 
complete list of systems with Vy < 6, supposing for the moment that 
we are not yet interested in the individual multiplicities. There are, 
however, minor issues, which make for an uncertainty of perhaps 
+10 in our list of 4559 systems. The main one is deciding whether 
two or more stars, each fainter than Vy = 6, are a real system or 
not. We believe that there is no answer to this that everyone would 
accept; and the main reason for this is the long-range character 
of the gravitational force, something that we cannot vary or work 
around. Nevertheless, we have tested a number of algorithms in 
which brightness is added up for all Hipparcos targets that are closer 
together than some limit – either an angular limit, or a linear limit 
involving parallax. With appropriate choice of limit, we can pick 
up clusters as spread out as the Hyades, but we have concentrated 
on smaller limits, in particular 180 arcsec. This produces 130 pairs 
or higher multiples, in addition to many more that qualify only as 
single targets in HIP but as multiples in the HDMC catalogue. The 
great majority of the 130 are well-known ‘systems’ whose reality 
has been investigated over decades or centuries, and where there is 
rather little doubt of the reality or otherwise. But we always come 
across a few near the margin, wherever and however we might try 
to define the margin. The smallest angular separation that we accept 
as optical (among pairs of nearly equal brightness) is ~27 arcsec, 
for HIP 35210/35213 (HR 2764) and also for HIP 79043/79045 
(HR 6008/6009), and the two largest that we accept provisionally 
as real аге 7860 arcsec for а Cen (HR 5460) and ~7030 arcsec, for a 
PsA (HR 8728). The linear separation in о PsA is slightly more than 
0.25 pc, and we can wonder whether such a system can survive even 
for one full rotation. The age of œ PsA is estimated as 200 + 100 Myr 
by Di Folco et al. (2004). If the two main bodies have been moving 
apart steadily over that interval, i.e. not orbiting but escaping, they 
have been doing so at the very low velocity of 0.0025 km s^!, much 
less than the actual escape velocity. On the other hand, if they are 
orbiting, they have survived rather surprisingly for 2,20 orbits of 
~6 Myr each. But despite the uncertainty of such cases, we feel that 
only about 10 in (or not in) our entire sample are rather marginal 
either way. 

Detection efficiency, such as the difficulty of observing long- 
period spectroscopic binaries and short-period visual binaries, will 
make the distributions of multiplicity, and of periods, differ from 
the true distribution. One of us (PPE) proposes to investigate this 
further using a Monte Carlo procedure that allows one to (1) con- 
struct a magnitude-limited sample of systems with multiplicity one 
to eight according to some hypothetical distribution of masses, ages, 
multiplicity and periods, and (ii) to ‘theoretically observe’ the mem- 
bers of this sample, making estimates of the efficiency of various 
observational procedures. This can map the original multiplicity 
into at least a lower limit to the observed multiplicity. It is difficult 
to believe that observations to date can constrain the number of T 
dwarf companions that an O star might have, but some less extreme 
pairings may be capable of being ruled in (or out). 

Selection effects often interact with detection efficiency: if a given 
technique is known to work well for certain types of star, those types 
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of star are more likely to be investigated with that technique. Thus, 
excellent radial-velocity orbits of quite long-period G/K binaries, 
say 1-30 yr, are known, but very few in this period range are known 
for O or B stars because their broad lines preclude measurement 
to the necessary accuracy. Nevertheless, there may be the same 
proportion of 1—30 yr binaries in both subsamples. 

While the limits of various observing techniques in detecting bi- 
nary companions can be modelled reasonably well, the extent of 
their application to our sample remains unexplored. For example, 
radial velocities of almost all bright stars have been measured sev- 
eral times, but the accuracy and time coverage vary to such an extent 
that we cannot apply uniform criteria to the sample as a whole. Many 
subsystems which are detectable spectroscopically still remain un- 
detected. Similarly, only a fraction of bright stars has been observed 
interferometrically. 

We note in Table 5 that the period distribution is bimodal at 
early spectral types, but unimodal, with a very broad maximum, 
at later types. That might be because the intermediate periods for 
O stars, 0.1—100 yr, are too close to recognize visually (although 
some are resolved), and too wide to recognize spectroscopically 
(although some also have been recognized). A preliminary version 
ofthe Monte Carlo code (Eggleton, Kisseleva-Eggleton & Dearborn 
2007) suggested that the apparent bimodality of O-star periods can 
be explained this way; we should bear in mind the result (Table 6) 
that O stars appear to be more highly multiple than later stars, so 
that it can be not improbable for an O-star system to have at least 
three bodies, one pair of which might have a period in the ‘missing’ 
range. 

There are indications in Table 4, or equivalently Fig. 1, that many 
multiples remain to be detected. On the one hand, the 4-mag-limited 
sample of Table 4 shows frequency N, of multiplicity n dropping 
substantially less rapidly (roughly, №, ~ 2.37") than the 6-mag- 
limited sample (roughly, №, ~ 3.47"). We would expect at most 
a small effect here. The brighter sample is on the whole nearer, 
so that visual multiples are more readily recognized, and biased to 
slightly more massive stars, which tend to be more highly multiple. 
But preliminary attempts at Monte Carlo modelling (Eggleton et al. 
2007) suggest that this can only account for perhaps a third of the 
change in average multiplicity as seen between the two samples. We 
suggest that most of the difference is due to the fact that over the last 
three centuries the brightest stars have been studied most carefully. 


n 


log N 


Figure 1. Log frequency as a function of multiplicity, from the first two 
lines of Table 1. Least-square fits are plotted, with slopes corresponding to 
—3.37 for the larger sample and —2.30 for the smaller sample. 
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Interestingly, this may now be the opposite of the truth. Modern 
detectors can so be sensitive that observers avoid the brightest stars. 

The comparison of the northern and southern half-samples in 
Table 4 tells a possibly similar story. The fact that the northern 
sample has a modestly higher average multiplicity than the southern 
sample may be because the number of telescopes, and of observers, 
in the northern hemisphere has always been substantially larger. 
The difference between the northern and southern samples is not 
significant at the 99 per cent level, but it is significant (C. A. Booker 
2008, p.c.) at the 95 per cent level, according to a x? test. The 
preliminary Monte Carlo model of Eggleton et al. (2007) suggested 
that the true average multiplicity would have to be above 2.0 if the 
observed average multiplicity for the complete sample (top line of 
Table 2) is to be in excess of 1.5. 

In Table 5, we do not discriminate between inner and outer pe- 
riods, although it would not be difficult in principle to do so. Such 
discrimination would reveal that close pairs prefer to be inner com- 
ponents of higher multiplicity systems rather than pure n — 2 bi- 
naries. The present catalogue, as well as the MSC, contains rich 
information on the statistics of periods and mass ratios at different 
hierarchical levels. 


6 DISCUSSION 


We see three main areas where statistics of the multiplicity and 
period distributions for a complete magnitude-limited sample is 
potentially useful. First, they can constrain modelling of the star 
formation process. Star formation lacks any clear initial conditions, 
but at least it has some reasonably clear ‘final conditions’, and 
models can be tested against their ability to produce the kind of 
systems that are actually seen, in the right proportions. Secondly, 
the model can serve as the initial condition for subsequent nuclear 
evolution; and as we amplify below there is also the possibility of 
dynamical evolution within triples and higher multiples that can 
influence their nuclear evolution. Thirdly, we believe that such a 
model can also be helpful in computing the evolution of star clusters 
under the joint effects of nuclear evolution and N-body gravitational 
dynamics. Although the dense environment of those star-forming 
regions that generate clusters, i.e. that make a transition from a gas- 
dominated entity to a star-dominated entity of much the same size, 
is likely to prevent the existence of the widest multiples that our 
sample contains, it is a reasonable first hypothesis that the shorter 
period multiples will not be very different. If, in fact, an N-body 
calculation was started with statistical multiplicity as in our sample, 
presumably the wider systems would be rapidly broken up; but 
the closer systems could persist, and the presence of rather close 
primordial triples and quadruples would influence significantly the 
evolution of the cluster. 

Regarding the second area, Eggleton & Kiseleva (1996) have 
enumerated a number of ways in which the existence of primordial 
triples can allow evolutionary channels that are different from what 
can be obtained from only binaries. We mention two here. First, 
triple stars in which the two orbits are misaligned can be subject 
to the dynamical effect of Kozai cycles (Kozai 1962), and these in 
turn can allow tidal friction to become important in the course of 
105—10? yr and cause the inner orbit to become smaller (Mazeh & 
Shaham 1979; Kiseleva, Eggleton & Mikkola 1998; Eggleton & 
Kisseleva-Eggleton 2006). We call this process Kozai Cycles with 
Tidal Friction (KCTF). Tokovinin et al. (2006) have noted, using 
a maximum-likelihood method, that in a sample of spectroscopic 
binaries with periods «2.9 d, 96 + 7 per cent have third bodies, 
as compared with 34 + 6 per cent in the period range 13-30 d. 


Fabrycky & Tremaine (2007) have computed the effect of KCTF 
on a Monte Carlo ensemble of triples, and shown that indeed it can 
produce an accumulation at short inner periods. Pribulla & Rucinski 
(2006) have noted that as many as 42 per cent of contact binaries 
appear to be in triples (and arguably 59 per cent in a more thoroughly 
examined subsample), and it could be that KCTF has contributed to 
this; although we probably need the additional influence of magnetic 
braking, also with tidal friction (MBTF), to drive fairly close low- 
mass binaries generated by KCTF to contact on a time-scale of 
109. In our primary sample, line 1 of Table 4, there are 95 binaries 
with P < 3d, and 60 of these are in systems with n > 3; this is a 
much higher proportion than for longer periods. 

The second effect is the production of ‘anomalous binaries’. In 
short-period binaries, we can expect that a merger of the two com- 
ponents is a fairly common event. Case A systems can evolve con- 
servatively only if the initial mass ratio is fairly mild (1 > q 2, 0.6; 
Nelson & Eggleton 2001), and if q is not this mild then a merger 
seems quite a likely event. It would be hard to determine that a 
particular currently single star is a merged remnant of a former 
binary, although some Be stars that are apparently single might be 
such remnants. But within a primordial triple system, it is possible 
that such a merged remnant would be identifiable, because the wide 
binary that remains after the merger of the close subbinary would be 
expected, in at least some circumstances, to show an anomaly where 
the two components appear to be of different ages. R. E. M. Griffin 
(to be published) has found a number of such apparently anomalous 
systems, of which y Per (HR 915) is an example. Although the 
mass ratio, obtained from careful deconvolution of the two spectra 
(G8II-III + A2IV; 53304, e = .79), is 1.54 (in the sense Ma/MA), 
the A component seems surprisingly large and luminous compared 
with what it should be on the ZAMS; and it ought to be very close 
to the ZAMS if itis coeval with the G component. A possible expla- 
nation is that the giant is the merged product of a former subbinary 
with a mass ratio of ~0.5, since this could allow the more mas- 
sive two of the original three components to evolve at roughly the 
same rate (Eggleton & Kiseleva 1996). We hope to test shortly the 
possibility that the appropriate primordial triple parameters, from 
our Monte Carlo model, will give an acceptable number of poten- 
tial progenitors. Alcock et al. (1999) and Evans et al. (2005) have 
noted that a similar process might lead to Cepheid binaries of an 
anomalous character, such as may be required to reconcile observed 
Cepheids with the theoretical models of the Cepheid pulsation 
phenomenon. 
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